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HTL quasiparticle picture of the thermodynamics of QCD* 

A. Rebhan a 

a Institut fiir Theoretische Physik, Technische Universitat Wien 
Wiedner Hauptstr. 8-10, A-1040 Vienna, Austria 

Starting from a nonperturbative expression for entropy and density obtained from ^-derivable two-loop ap- 
proximations to the thermodynamic potential, a quasiparticle model for the thermodynamics of QCD can be 
developed which incorporates the physics of hard thermal loops and leads to a reorganization of the otherwise 
ill-behaved thermal perturbation theory through order a^ 2 . Some details of this reorganization are discussed 
and the differences to simpler quasiparticle models highlighted. A comparison with available lattice data shows 
remarkable agreement down to temperatures of ~ 2.5T C . 
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1. Entropy and density in ^-derivable ap- 
proximations 

In the so-called ^-derivable approximations |Q 
obtained by truncating the skeleton functional 
($) in the Luttinger-Ward representation of the 
thermodynamic potential the expression for 
the entropy density takes a remarkable form: for 
example in a scalar theory gip 3 
written as 

d k dn^E^) 



g 2 <p it can be 



S = - 



(2tt) 4 dT 



■Im logD _1 (o;,fc) 



d A k dn(ui) 
(2tt) 4 dT 



Imll(cj, k) Re D(u, k) +S' (1) 



where D and II are self-consistent propagator and 
self-energy, resp., and one finds 

S' = (9(3-loop) = 0{g A ). (2) 

The fact that S' = in 2-loop <&-derivable ap- 
proximations, which has been observed first in 
Fermi- liquid theory by Riedel ||, turns out to 
hold rather generally || |^|. In particular, the 2- 
loop contribution of Dirac fermions to the entropy 
density reads 

d A k dnpD(uj) 



S 



^Mtrjlmlog^- 1 ) 



Im(7oE)Re(5 7 o)k (3) 



For nonvanishing chemical potentials //$ , this can 
be extended to the fermion number densities 
Mi, which are obtained by replacing the explicit 
derivative d/dT in (H) by d/d/Ai, with 



S' f = 0(3-loop) = M[ 



(4) 



Explicit interactions at 2-loop order are thus 
completely encoded in the dressed (quasiparticle) 
propagators D, S, . . ., which makes entropy and 
density a preferred starting point for a quasipar- 
ticle description of thermodynamic quantities — 
in fact any thermodynamic quantity, since apart 
from a single integration constant (which in QCD 
corresponds to the inherently nonperturbative 
bag constant) the grand canonical potential can 
be reconstructed from entropy and density. 

If a quasiparticle picture in the sense of Lan- 
dau is applicable, that is, if a large fraction of 
the possibly strong interactions can be accounted 
for by 'propagator renormalization', the above 
dressed 2-loop expressions should provide a rea- 
sonable approximation to thermodynamic quan- 
tities and the residual quasiparticle interactions 
~ g A described by S' or M' should be compara- 
tively weak. 

Technically, an important point turns out to be 
that the above expressions for entropy and den- 
sity are manifestly ultraviolet finite, as all inte- 
grals involve dn/dT or dn/d^i which go to zero 
exponentially for both positive and negative fre- 
quencies. In relativistic field theories, this is a 
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great advantage, because it allows one to use 
these functionals in a nonperturbative manner 
once finite approximations for the propagators 
and self-energies have been obtained, whereas 
usually the renormalization programme requires 
systematic expansions (and truncations) in pow- 
ers of the coupling constant .0 

At finite temperature, truncated perturbative 
series turn out to have extremely poor appar- 
ent convergence behaviour and this problem 
arises with the appearance of odd powers in the 
coupling such as g 3 (the so-called plasmon ef- 
fect^]). In QCD this plasmon effect spoils ap- 
parent convergence up to extremely high tem- 
peratures ^ 10 5 T C . However, this breakdown 
of perturbation theory does not necessarily have 
anything to do with the nonabelian character of 
QCD, as a similar phenomenon occurs in virtu- 
ally any field theory for coupling constants where 
T = perturbation theory still looks applicable. 

On the other hand, in expressions such as ([!]) 
and @ the coupling constants are all contained 
in dressed propagators and self-energies, which 
themselves are buried in logarithmic or fractional 
expressions. If there are contributions which 
when expanded out in powers of the coupling have 
uncomfortably large coefficients, they do not nec- 
essarily lead to large corrections as they are kept 
together with higher-order terms (which usually 
would be discarded) such that final results be- 
have more like logarithms or fractions rather than 
polynomials. 

Starting from the 2-loop ^-derivable entropy 
(density), one should therefore be in a much more 
favourable position to investigate whether and 
to what extent a quasiparticle picture such as 
that underlying hard-thermal- loop (HTL) pertur- 
bation theory jl2| is applicable. 

In the following I shall describe the results ob- 
tained by employing HTL propagators and self- 
energies and certain corrections thereof in the 
above nonperturbative expressions for entropy 



1 In the recent work of Ref. the renormalizability of 
<J?-derivable approximations has been postulated, which is 
however nossih lv i n contradiction with the explicit calcu- 
lations of Ref. [0. 

2 Caused in fact by the appearance of screening (Debye) 
masses as opposed to the dynamical plasmon mass. 



(and density). 

This will only provide an approximation to 
the self-consistent 1-loop gap equation that per- 
tains to the ^-derivable 2-loop approximation. A 
complete solution of the self-consistent gap equa- 
tion is presumably prohibitively complicated, but 
apart from that, not really desirable: While eqs. 
(|l|) and (||) are UV finite, the associated 1-loop 
gap equations are not. Moreover, in gauge the- 
ories these equations are gauge dependent, as 
<I>-derivable approximations do not respect lo- 
cal gauge invariance. The approximately self- 
consistent solutions that will be discussed below 
evade both problems by providing finite, gauge 
invariant and gauge fixing independent approxi- 
mations, and they do so by dropping terms that 
are anyway beyond the accuracy of any 2-loop 
approximation, i.e., of order <? 4 or higher. 



2. Approximate propagators 

In thermal perturbation theory one has to dis- 
tinguish hard and soft momentum scales. At high 
temperature (in QCD: sufficiently above the de- 
confinement transition), the dominant contribu- 
tion to the grand canonical potential comes from 
elementary quanta with energies and momenta 
~ T, which will be called "hard" . Collective phe- 
nomena arise at the scale gT, which will be called 
"soft". Even when g ~ 1, as will be the case of 
interest in QCD, this distinction may still make 
sense, if the scale of hard physics is actually 2itT 
(the spacing of the Matsubara frequencies)rather 
than just T, as frequently argued jlql . 

At soft momentum scales, the leading terms 
in the self-energies are given by the HTL self- 
energies, which for gauge bosons and (massless) 
fermions contain two separate structure functions 
each and read [Hi 



n L (w, k) = m 2 D 



^ US j us + k 

2k us — k 



E±(us,k) 



U T (u,k) -- 
M 2 



US 



us =F k ^ us + k s 
2k ^ us — k 



(5) 
(6) 
(7) 
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with 



2tt 2 ' 



(8) 



(9) 



The gauge boson propagator involves two dif- 
ferent propagators, a spatially transverse one, 

D t (uj, k) = [-uj 2 + k 2 + ft x (w, fe)]- 1 , (10) 

and one that only appears in the thermal 
medium, 



D L (u,k) = -[k 2 + IL L (u,k)]- 1 . 



(11) 



Similarly, the fermion propagator has a "normal" 
branch S+ cx [uj — (k + and one that is 

a purely collective effect, SL cx [uj + (k + 
where the label (— ) refers to the fact that chirality 
and helicity have opposite signs. 

These HTL propagators have poles at real fre- 
quencies uj for all k, corresponding to (undamped) 
quasiparticles with momentum-dependent ther- 
mal masses, with the additional branches disap- 
pearing from the spectrum for k/gT 3> 1 because 
of exponentially vanishing residues. The normal 
branches on the other hand, approach constant 
asymptotic masses m 2 ^ = ^rh 2 ^ and = 2Af 2 , 
resp., in this limit. 

At hard momenta, the above HTL expressions 
are actually invalid, except for \uj 2 — k 2 \ <C T 2 . 
This is however the region where all the poles of 
the HTL propagators lie, so that the latter still 
give the correct dispersion law of quasiparticles 
at leading order, with next-to-leading order cor- 
rections calculable within standard HTL pertur- 
bation theory. 

2.1. Leading-order terms in entropy and 
density 

Let us now inspect how the leading-order in- 
teraction terms of the grand canonical potential 
arise in the self-consistent entropy (and density) 
functional (|l]), @, which at two-loop order are 
determined by interaction-free quasiparticle prop- 
agators. 

Since H cc g 2 , these must be reproduced by 
linearizing entropy and density in n. In the pure 



glue case, the relevant terms in the entropy read 

T (2) 



d 4 k &n r n T 

— — < — Im 



(2tt) 4 dT 



■Imll^Re 



uj 2 — k 2 
1 



= 2N„ 



d k dn 



(2tt) 4 dT 



^— Reni 2) Im 



1 



uj 2 — k 2 ' 

-7Tc(w)5(^ 2 -fc 2 ) 



(12) 



where only the g 2 part of the transverse polariza- 

(2) 

tion tensor, T1 T is needed. The latter is seen to 
contribute only its light-cone value 



Tl {2 \uj 2 = k 2 ) = g 2 NT 2 /6 = m l, 



(13) 




yielding S 2 = —Ngin^T/Q with N g = N 2 - 1. 

Analogous results are obtained for the 
fcrmionic contributions, and also in the case of 
fermion densities. The order-g 2 terms are simply 
and generally given by the (leading-order) ther- 
mal masses of hard particles according to 

Kb f 

F 

where the sums run over all the bosonic (B) 
and fcrmionic (F) degrees of freedom (explicitly 
counting spin degrees of freedom). 

2.2. Plasmon terms 

In conventionally resummed perturbation the- 
ory, the so-called plasmon term cx g 3 in the ther- 
modynamic potential arises from the appearance 
of the Debye mass in the electrostatic propaga- 
tor. Only the Matsubara zero modes are able to 
contribute an odd power of the Debye mass (and 
thus an odd power in the coupling) through ring 
resummation according to 



P 3 = -NgT 



d 3 k 
(2^)3 



log(l + ^)-^ 



m 2 D 



= N„ 



12tt 



(16) 



In the self-consistent entropy expression, this 
plasmon term arises rather differently. It no 
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longer suffices to resum the Debye mass for zero 
modes. Instead, there are soft as well as hard 
contributions to order g 3 . 

At soft momenta, these contributions come 
from HTL propagators and self-energies in the 
form (for pure glue and with implicit sums over 
colour and polarization states) 



^soft 



d A k 1 



2*) 



-{ Im 



io g (i + D n) - nD 



-ImIIRe(£> - D Q 



dP 3 
dT 



with 



AS 3 = N„ 



AS, 



d 4 k 1 



(2tt) 4 uj 



21m 



n T Re (b T 



- D 



(17) 



(0) 



Imn L Rc [D L - D 



>(o) 



AS^ + AS^ 



(18) 



We found numerically that AS3 — by cancel- 
lations in more than 8 significant digits, without 
being able to gain more fundamental insight into 



the peculiar sum rule AS^ ' = —AS 



:(3) 



which 



emerges only after carrying out both, the fre- 
quency and the momentum integrations in fll8|). 
In the pure-glue case we thus have 



5^ 



dP 3 
dT 



-1*. 



(19) 



At hard momenta, the only possibility for con- 
tributions ~ g 3 is through NLO corrections to the 
spectral properties of the hard excitations. In- 
deed, one can show that the remaining three 
quarters of the plasmon effect are provided by 



hard 



-N, 



d 3 k 1 dn(k) 



(2tt) 3 k dT 



ReSn T (tu = k). (20) 



HeSHxiuj = k) is a momentum-dependent cor- 
rection ~ g 3 T 2 to the asymptotic thermal mass 
of transverse gluons, and it is calculable from 
standard HTL perturbation theory without be- 
ing afflicted by the IR sensitivity that arises in 
the imaginary part jl5| . 

Clearly, the self-consistent entropy involves a 
"massive" reorganization of thermal perturbation 



theory Instead of coming exclusively from a mod- 
ification of the static propagators, the plasmon 
effect is now spread over the entire spectral prop- 
erties of quasiparticles as given by HTL propaga- 
tors at soft momenta and HTL masses plus NLO 
corrections at hard momenta. While this seems 
to be just an extravagant complication from the 
point of view of perturbation theory, with regard 
to the self-consistent entropy it is gratifying that 
the plasmon effect is encoded in the spectral data 
of (HTL) quasiparticles in a more uniform man- 
ner, with the bulk of the plasmon effect coming in 
fact from the dominant hard degrees of freedom. 

All this can be extended in a straightforward 
manner to Nf 7^ 0. In the case of finite chemical 
potential but high enough temperature T ^> rhu, 
one finds that more than 1/4 of S3 is coming from 
the soft momentum region, whereas all of the 
plasmon term in the density, A/3, is due to NLO 
corrections of the asymptotic thermal fermion 
masses. 

3. Nonperturbative usage 

As mentioned above, conventional thermal per- 
turbation theory exhibits very poor apparent con- 
vergence as soon as collective phenomena such 
as the plasmon effect are included, which involve 
odd powers in the coupling. The self-consistent 
entropy and density functionals at two-loop or- 
der also contain the plasmon effect, but together 
with higher-order terms that ordinarily would 
have been discarded by a strictly perturbative 
expansion in g because otherwise the renormal- 
ization programme could not have been carried 
out. By contrast, the two-loop entropy and den- 
sity functional is UV finite and can be evalu- 
ated in a nonperturbative manner using the high- 
temperature/density approximation of the gluon 
and quark propagators, which are finite at lead- 
ing and next-to-leading order in HTL perturba- 
tion theory. 

3.1. HTL/HDL approximation 

As a first approximation let us consider the 
two-loop entropy and density functionals evalu- 
ated completely using HTL (HDL) propagators. 
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3.1.1. Entropy 

In the case of the entropy, this takes care of all 
contributions of order g 2 , but only part (1/4) of 
the plasmon term ~ g 3 . However, it also contains 
infinitely many higher-order terms which despite 
being incomplete may help to get rid of the patho- 
logical behaviour of the perturbation series trun- 
cated at low orders in g. Among such higher- 
order contributions is for instance a <? 4 contribu- 
tion reading (for pure glue) 

^—"•Sf{^ + 1M ") (21 » 

involving a g 4 log(l/<7)-term, whose coefficient in 
pure-glue QCD is 1/12 of that of the complete 
perturbative result.]] By contrast, simple mas- 
sive quasiparticle models such as those used in 
Refs. |l6| , |l7| do not have a g 4 log (1/g) -term in 
the entropy at all. 

In Fig. [j], the numerical result for Shtl/Ssb 
in the case of pure glue is given as function of 
tHd/T, which is the only independent parame- 
ter. The HTL result is given by the full line 
and is found to be a monotonicly decreasing func- 
tion of tHd/T. If this were expanded in powers 
of ijid/T and truncated beyond (ttid/T) 3 ~ g 3 
(dashed line in Fig. |l|), this property would have 
been lost at rhry / (2ttT) w 1/3, where one might 
still expect a sufficiently clear separation of hard 
and soft scales which is a prerequisite of the HTL 
approximation. 

The numerical result for Shtl is in fact very 
close to a simple massive quasiparticle model with 
entropy 2N g SsB(jn 00 ), represented by the dot- 
ted line in Fig. |l|. Remarkably one has Shtl < 
2NgSsB (toqo) even though the latter contains 
30% less of the plasmon term ~ g 3 , which would 
be expected to work in the opposite direction. 
This is further emphasized by comparison with 
its perturbative approximation given by the dash- 
dotted line. 

The rather small difference to a simple mas- 
sive quasiparticle model is further inspected in 
Fig. |[ where it is shown that there are large can- 

3 The correct coefficient will be restored by 
0(cj 4 log(l/</)T 2 ) corrections to m 2 ^, whereas the 
constant under the logarithm also receives three-loop 
contributions. 
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Figure 1 . HTL entropy per gluonic degree of free- 
dom normalized to its Stefan-Boltzmann value as 
a function of the Debye mass ttid(T, ll)/T. The 
full line gives the complete numerical result; the 
dashed line corresponds to the perturbative re- 
sult to order (m D /T) 3 ~ g 3 . The dotted line 
gives the entropy for scalar degrees of freedom 
with momentum-independent mass m = TOqo = 
tHd/ a/2; its perturbative approximant is given by 
the dash-dotted line. 



cellations between the longitudinal and transverse 
contributions to the HTL entropy. 

In Fig. H the HTL entropy of 1 quark de- 
gree of freedom at zero chemical potential is dis- 
played as a function of the fermionic plasma fre- 
quency M/T. Like in the pure-glue case, this 
contains also an (incomplete) <? 4 log(l/g) contri- 
butions that is not present in simpler quasiparti- 
cle models: 

Sf H TL = NN f^f I + °' 22 ' ' ■) (22) 

Sf.HTL does not contain any contributions to the 
plasmon term ~ g 3 , which entirely come from 
NLO corrections to M^. There is therefore no 
big deviation of the full numerical result (solid 
line) from the perturbative one truncated beyond 
g 2 (dashed line). 

Compared to the entropy of a simple massive 
fermionic quasiparticle model S^(Moo) (dotted 
line in Fig. ||), one finds extremely good numeri- 
cal agreement, which again takes place only after 



G 



[Shtl - 2N g S ( mao )]/\2N g S (m 




Figure 2. Relative deviation of the HTL entropy 
from that of a gas of massive bosons with (con- 
stant) mass moo (full line) and their composition 
from transverse (T) and longitudinal (L) quasi- 
particle (QP) and Landau damping (LD) contri- 
butions. 
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Figure 3. HTL entropy per quark degree of free- 
dom at fi — normalized to its free value (solid 
line), the corresponding perturbative order-g 2 re- 
sult (dashed line), and the entropy of a quark 
with constant mass (dotted line). 



adding up all the quasiparticlc ((+) and (— )) and 
Landau-damping contributions, however. 

3.1.2. Quark density 

In the case of zero temperature but high chemi- 
cal potential, the quantity of interest is the quark 
density. This does not contain any plasmon term 
~ g 3 , but rather ~ 3 4 log(l/.g). The HDL ap- 
proximation does contain some though not all of 
this term: 



M 4 



log 4- + 0.35 . 
M 



(23) 



Ordcr-g 2 log g contributions to the asymptotic 
masses of the quark and gluon quasiparticlcs, still 
within the framework of the 2-loop ^-derivable 
approximation, are responsible for the remaining 
contribution to the coefficient of the g 4 log (/-term, 
while the coefficient under the logarithm also re- 
ceives 3-loop contributions. 

In Fig. H the numerical result for Mh dl is given 
as a function of M / fi (full line) and compared to 
that of a simple quasiparticlc model 



T=0 



1 

3T 2 



(fi 2 -M^e^-M 0O )(24) 



as well as to a perturbative approximation trun- 
cated beyond (M//J,) 2 ~ g 2 . Remarkably, the 
full HDL result drops to zero at almost the same 
value as the simple quasiparticle model. However, 
the former becomes negative thereafter, showing 
that the HDL approximation is breaking down at 
Moo/ A* > 1 at the latest. 

3.2. Next-to-leading approximations 

The plasmon term ~ g 3 at high temperatures 
T rfi£) becomes complete only after inclusion 
of the next-to-leading correction to the asymp- 
totic thermal masses and M^. These are 
determined in standard HTL perturbation theory 
through 

Smlcik) = RedTl T {uj = k) 

-A 



Re( 



(25) 



where thick dashed and wiggly lines with a blob 
represent HTL propagators for longitudinal and 
transverse polarizations, respectively. Similarly, 

(26) 



Re( 



The explicit proof that these contributions in- 
deed restore the correct plasmon term is given in 
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Figure 4. HDL quark density at T = normalized 
to its free value (solid line), the corresponding 
perturbative order-g 2 result (dashed line), and 
the free quark density of a quark with constant 
mass Mqo (dotted line). 



Rcf. §. 

These corrections to the asymptotic thermal 
masses are, in contrast to the latter, nontrivial 
functions of the momentum, which can be evalu- 
ated only numerically. However, as far as the gen- 
eration of the plasmon term is concerned, these 
functions contribute in the averaged form 



5ml 



j dkkn' BE (k)RedU T {uj = k) 
J dkk n' BE (k) 



(cp. eq. |l4| ) and similarly 

2 _ / dkkn' Frt (k)Rc2kdY,+ {uj = k) 



5Mi 



j dkk np D (fc) 



(27) 



(28) 



These averaged asymptotic thermal masses turn 
out to be given by the remarkably simple expres- 
sions 



<5m_ = 



1 



NTrhi 



6M °° = -^C f Tm D , 



(29) 

C f = N g /(2N). (30) 



These corrections only pertain to the hard 
excitations; corrections to the various thermal 
masses of soft excitations are known to differ sub- 
stantially from ( p9| ) or d3(i|). For instance, the 
gluonic plasma frequency at k = reads JL8[ 



brr? pl /rhpi ~ — 0.18y/Ng, which is only about a 
third of Jm^/m^; the NLO correction to the 
nonabelian Debye mass on the other hand is even 
positive for small coupling and moreover logarith- 
mically enhanced (T^j 8m 2 D /rh\) = +\/3N /(2n) x 
Slog(l/s). 

For an estimate of the effects of a proper incor- 
poration of the next-to-leading order corrections 
we have therefore proposed to include the latter 
only for hard excitations and to define our next- 
to-leading approximation (for gluons) through 



Snla — Shtl + Sfn'. 

soft 30 hard 



(31) 



where m 2 ^ includes (|2^). To separate soft {k ~ 
rho) and hard (k ~ 2nT) momentum scales, we 
introduce the intermediate scale A = ^/2irT , m£,c^ 
and consider a variation of ca = \ ■ ■ ■ 2 as part of 
our theoretical uncertainty. 

Another crucial issue concerns the definition 
of the corrected asymptotic mass moo. For the 
range of coupling constants of interest (g > 1), 
the correction l^m^l is greater than the LO value 
ttt-oo , leading to tachyonic masses if included in a 
strictly perturbative manner. 

However, this problem is not at all specific to 
QCD. In the simple g 2 (p A model, one-loop re- 
summed perturbation theory gives 



m 2 =g 2 ni--g) 

7T 



(32) 



which also turns tachyonic for g > 1. On the 
other hand, the solution of the corresponding sim- 
ple scalar gap equation is a monotonic function in 
g, and it turns out that the first two terms in a 
(m/T)-expansion of this gap equation, 



g ^T 2 - -g 2 Tm 

IT 



(33) 



which is perturbatively equivalent to ( p2|) has a 
solution that is extremely close to that of the full 
gap equation (for MS renormalization scales p, ~ 
2ttT). 

In QCD, the (non-local) gap equations are way 
too complicated to be attacked directly. We in- 
stead consider perturbatively equivalent expres- 
sions for the corrected moo which are monotonic 
functions in g. Besides the solution to a quadratic 
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Figure 5. Comparison of the lattice data for 
the entropy of pure-glue SU(3) gauge theory of 
Ref. [M (gray band) with the range of Shtl 
(solid lines) and Smla (dash-dotted lines) for 
P = ttT... AttT and c A = 1/2 . . . 2. 




T/T c 

Figure 6. Comparison of the HTL entropy (solid 
lines) and the NLA estimate (dash-dotted lines) 
for Nf = 0,2,3 with the estimated continuum 
extrapolation of Nf = 2 lattice data of Ref. p2{ . 



equation analogous to (33) we have tried the sim- 
plest Pade approximant m 2 = g 2 T 2 /(l + ^g), 
which also gives a greatly improved approxima- 
tion to the solution of scalar gap equations. In 
QCD, our final results do not depend too much 
on whether we use the Pade approximant (||| or 
a quadratic gap equation |7). 

The main uncertainty rather comes from the 
choice of the renormalization scale which deter- 
mines the magnitude of the strong coupling con- 
stant when this is taken as determined by the 
renormalization group equation (2-loop in the fol- 
lowing) . 

In Fig. |, the numerical results for the HTL en- 
tropy and the NLA one are given as a function of 
T/T c with T c chosen as T c = l.UA^. The full 
lines show the range of results for Shtl when the 
renormalization scale ft is varied from ttT to AttT; 
the dash-dotted lines mark the corresponding re- 
sults for Snla with the additional variation of 
ca from 1/2 to 2. The dark-gray band are lattice 
data from Ref. |2(if] ; the more recent results from 
Ref. 1 21 1 are consistent with the former within 



error bars and centered around the upper bound- 
ary of the gray band for T w 3T C and somewhat 
flatter around 2T C . Evidently, there is very good 
agreement for T > 2.5T C . 



In Fig. |6[ Nf massless quarks are included and 
compared with a recent estimate [ p2[ of the con- 
tinuum limit of lattice results for Nf = 2 (gray 
band), but now with Shtl and S^la evaluated 
for the central choice of /2 = 27rT and ca = 1 
(with unchanged T c /Aj^). When Nf is in- 
creased, there are competing effects of larger ther- 
mal masses versus slower running of a s , which re- 
sult into a rather weak dependence of our results 
on Nf as a function of T/Aj^-j as it is in Fig. [| 

From the above results for the entropy density, 
one can recover the thermodynamic pressure by 
simple integration, 



P(T) - P{T ) = [ dT'S{T'). 
Jt„ 



(34) 



The integration constant P(Tq), however, is a 
strictly nonperturbative input. It cannot be fixed 
by requiring P(T = 0) = 0, as this is in the con- 
finement regime. It is also not sufficient to know 
that limT^oo P = -Pfree by asymptotic freedom. 
In fact, the undetermined integration constant in 
P(T)/Pf lee (T) when expressed as a function of 
a s {T) corresponds to a term B 



C exp 



1 



(35) 
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Figure 7. Comparison of our results for the pres- 
sure of QCD with Nf — 2 with the extrapolated 
lattice data of Ref . |E3| . 



which vanishes for a s — > with all derivatives 
and thus is not fixed by any order of perturbation 
theory. It is, in essence, the nonperturbative bag 
constant, which can be added on to standard per- 
turbative results, too. However, in P(T)/ Pf Tee (T) 
this term becomes rapidly unimportant as the 
temperature is increased, as it decays like T~ 4 . 

In Fig. ^ which shows the HTL and NLA pres- 
sure in the case of Nf = 2, we have fixed this 
integration constant by choosing P{T C ) = as 
T = T c is beyond the range of applicability of any 
form of perturbation theory anyway. As with the 
entropy, there is good agreement with lattice data 
for T > 2.5T C , and this agreement does not de- 
pend on the precise value of P(T C ), which enters 
only at the percent level at such temperatures. 

4. Discussion and outlook 

Starting from the self-consistent two-loop ap- 
proximations to entropy and density, we have 
been able to give a quasiparticle description of the 
thermodynamics of QCD that agrees well with 
available lattice data for T > 2.5T C , and which 
incorporate more of the perturbatively accessi- 
ble details of gluonic and fermionic quasiparti- 
cles than previously proposed quasiparticle mod- 
els |l6|,|l7|,|3| involving simply free massive scalars 
and fermions. In contrast to the latter, however, 



we refrained from introducing phenomcnological 
functions to fit the data, the only input being 
T c /Aj^-j as gleaned from the lattice. Put conser- 
vatively, we obtained a lower bound on T/T c for 
which HTL quasiparticles may be an appropri- 
ate description of hot QCD. This lower bound is 
however extremely low when compared with what 
is needed for conventional thermal perturbation 
theory to avoid an obvious breakdown, raising 
the hope that the latter just needs further resum- 
mations to become applicable at temperatures of 
physical interest in QCD. 

There exist in fact alternative proposals for 
such resummations. One particularly interesting 
approach is an extension of so-called screened per- 
turbation theory [^,^5|, where thermal masses 
are introduced as variational parameters, to non- 
abelian gauge theories, using the HTL effective 
action as gauge invariant thermal mass term p6| . 
It should be noted, however, that the latter is 
introduced as a mere technical device, whereas 
in the entropy-based approach the emphasis is 
on a description in terms of HTL quasiparticles 
and their perturbatively accessible modifications 
at NLO. Indeed, in the two-loop entropy it turned 
out that at hard momentum scales only the region 
close to the light-cone matters, where the HTL 
approximations remain accurate even at hard mo- 
menta. By contrast, in HTL-screened perturba- 
tion theory, the HTL's contribute throughout all 
of phase space. Moreover, the thus modified UV 
structure entails new UV problems (together with 
new scheme dependences). 

To avoid these UV problems, Peshier has 
proposed a somewhat contrived way of HTL re- 
summation directly on the thermodynamic pres- 
sure, which appears to come very close to what 
we have obtained in the HTL approximation in 
the entropy. Peshier's pressure is perturbatively 
equivalent to ours in the HTL approximation 
through at least order g 3 , though not obviously 
identical to it. 

Beyond the HTL approximation, it is not clear 
how the formulae of Ref. J27J could be extended. 
In the present HTL-based quasiparticle approach, 
the aim is to describe these quasiparticles in as 
full detail as HTL perturbation theory permits. 
In this respect we have so far taken only a first 
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step by considering averaged asymptotic thermal 
masses. We intend to study the full momentum 
dependence of the asymptotic thermal masses and 
how this enters into the thermodynamic quanti- 
ties in a future work. 

Another extension of the present work will be 
to consider in more detail the dependence on 
chemical potential, for example by integrating en- 
tropy and density for all T and [i in the deconfincd 
phase, similarly to what has been carried out in 
simple quasiparticle models P8fl . Using the NLA, 
this will include a complete plasmon term, and re- 
quires only one integration constant, P(To)| m= Oj 
that can be taken from the lattice. 

A first (infinitesimal) step towards finite \i can 
be studied on the lattice in the form of quark 
number susceptibilities dM /dfi\^ = o at high T 
p9| , p0| . Here conventional thermal perturbation 
theory is as ill-behaved as in the case of the pres- 
sure, whereas our approach can be applied in a 
straightforward manner, as will be presented in a 
forthcoming paper 
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